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Physics governing the locomotion of microor¬ 
ganisms and other microsystems is dominated 
by viscous damping. An effective swimming 
strategy involves the non-reciprocal and peri¬ 
odic deformations of the considered body. Here, 
we show that a magnetocapillary-driven self- 
assembly, composed of three soft ferromagnetic 
beads, is able to swim along a liquid-air inter¬ 
face when powered by an external magnetic field. 
More importantly, we demonstrate that trajec¬ 
tories can be fully controlled, opening ways to 
explore low Reynolds number swimming. This 
magnetocapillary system spontaneously forms by 
self-assembly, allowing miniaturization and other 
possible applications such as cargo transport or 
solvent flows. 

A Reynolds number much smaller than unity indicates 
that viscous forces dominate over inertial forces in a given 
flow. This is usually the case at the microscopic scale, 
which strongly impacts locomotion mechanisms of both 
biological and artificial microswimmers. Indeed, to swim 
by changing its shape, a microscopic body must break 
time-reversal symmetry m- Methods for producing 
artificial microswimmers are being actively researched 
Hi]. One major unanswered problem is the control of 
swimming direction. A fine control of swimming tra¬ 
jectories would be required for most practical applica¬ 
tions, such as manipulation and transport of small ele¬ 
ments [7] or micro-scale fluid flow generation [8] . Further¬ 
more, complex fabrication processes could limit potential 
micro-scale applications. Self-assembled systems present 
a clear advantage in this regard, as they require no direct 
manipulation of micro-components EHII]- Floating soft- 
ferromagnetic particles have been shown to produce self- 
assembled structures when exposed to magnetic fields 
m- Oscillating fields can deform these assemblies in 
a way that produces low Reynolds locomotion, although 
by which mechanism remains obscure m- In the present 
paper, we propose an explanation for the breaking of 
time-reversal symmetry in magnetocapillary swimmers, 
as well as demonstrate how swimming speed and direc¬ 
tion can be finely controlled in order to produce the de¬ 
sired trajectories. 

Three soft ferromagnetic spheres of size D are gen¬ 
tly placed along a water-air interface, at the center of 
a triaxial Helmholtz system. The weight of the parti¬ 
cles creates menisci that induce an attractive capillary 
interaction between neighboring beads [14], as shown in 
Figure [^a,b). Soft ferromagnetic particles are used, in 
which magnetic dipoles are induced along the direction of 
B with negligible hysteretic behavior. A magnetization 
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FIG. 1. Self-assembled magnetocapillary swimmer - 

(a) Sketch of the magnetocapillary system. Soft ferromag¬ 
netic particles are self-assembling at the water-air interface 
in a Petri dish. A vertical and constant magnetic field Bz 
is applied through the system. An oscillating horizontal field 
Bx excites the motion of the self-assembled structure, (b) 
A top view of the interface emphasizes the self-assembly of 
D = 500 fim beads in a vertical field : capillary attraction is 
counterbalanced by dipole-dipole repulsion. The deformation 
of the liquid-air interface is evidenced by placing an array of 
pixels underneath the Petri dish, (c) Five snapshots of the 
beads over one period T of the forcing oscillation. The traces 
of the initial positions are drawn in yellow for emphasizing 
the motion of the structure. One should notice the rotational 
oscillation of the structure during the period. 


cycle of the beads can be found at m- As a consequence, 
the breaking of time reversibility in our system does no 
come from magnetic hysteresis, as was proposed in recent 
theoretical models m- When the particles are in pres¬ 
ence of a magnetic field B = P BzCz^ dipole-dipole 
interaction competes with capillary attraction. A verti¬ 
cal magnetic field {Bx = 0) leads to the formation of a 
regular triangle, which size can be tuned by Bz IHIIT], 
as shown in the supplementary video 1. In the follow¬ 
ing, the vertical field is kept constant at = 30 G such 
that the bead center-to-center distance is close to two 
bead diameters 2D. When a constant horizontal field is 
added, symmetry breaking occurs in the system which 
then forms an isosceles. An oscillating horizontal field 
Bx = Px sin(27r/t) thus periodically deforms the trian¬ 
gular structure. Amplitudes ^x larger than Bzl2 cannot 
be reached because the system collapses, leading to hys¬ 
teretic contacts mi between the beads. 

In Figure [^c), one observes that small deformations 
of the triangle are accompanied by large rotational mo¬ 
tions of the whole structure. Those repeated deforma¬ 
tions/rotations come from the competition between mag- 
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FIG. 2. Selection of swimming modes -Probability Dis¬ 
tribution Function (PDF) of normalized speeds obtained in 
similar experimental conditions : I3x = 7.5 G and / = 0.5 Hz, 
but (a) without an offset, (b) with a constant offset Bxo = 
0.75 G. The plots are obtained from 55 independent realiza¬ 
tions. 


netic, capillary and hydrodynamic interactions m- A 
significant motion of the bead triplet can be seen in Fig¬ 
ure or in the supplementary video 2. Swimming has 
been found for frequencies / in between 0.1 Hz and 3 Hz. 
In the following paper, a frequency of 0.5 Hz has been 
chosen as a compromise between low Reynolds number 
and practical constraints linked to experiment running 
time. The relevant Reynolds number in the system is that 
of the individual particles, which in this case has typical 
values around Re ^ 10“^. It is possible, however, to fur¬ 
ther lower frequency / and/or field amplitude /3x in order 
to emphasize that low Reynolds locomotion takes place, 
reaching typical values in the range Re = [l0“^, 10“^]. 

Results 

Figure [^presents the Probability Distribution Function 
(PDF) of the normalized speeds for several swimmers 
using the same set of parameters : (a) without any offset 
and (b) with a small offset = 0.75 G. The swimming 
speed is found to be broadly distributed in the first case 
while it exhibits a narrow peak in the second case. In 
fact, various pulsation modes can be observed when the 
horizontal field oscillates. The offset Bxo helps the triplet 
to select a unique oscillation mode. The horizontal field 
becomes B^ = [^^co + sin(27r/t)]. One should remark 

that the offset enhances the deformation of the triangle 
but is kept as low as possible here in order to avoid the 
collapse of the structure. 

Figure [^a) shows the resulting speed v of the assem¬ 
bly normalized by the bead diameter D per forcing pe¬ 
riod T = 1//, as a function of the amplitude of the field 
oscillations /3x. For small oscillations (regime I), no sig¬ 
nificant motion is observed. When the amplitude /3x of 
the field oscillations becomes much larger than the offset 
Bxo (regime II), the speed increases and saturates at high 
field values. The dimensionless speed reaches about 0.3 
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FIG. 3. Swimming speeds - (a) Average speed of the self- 
assembled system as a function of the amplitude of oscilla¬ 
tions, when the offset is Bxo = 0.75 G and / = 0.5 Hz. The 
velocity is seen to increase and saturates around 0.3 DjT. 
The curve is a guide for the eye. (b) The maximum standard 
deviation S of the angles from the regular triangle as a func¬ 
tion of the amplitude /3x. The curve is a guide for the eye. 
Please note the large differences obtained near the collapse 
{/3x ~ 15 G). Three regimes (I, H and HI) are indicated and 
are discussed in the main text. 


D/T, which can be considered an efficient locomotion 
speed in the Stokes regime UHlEo]. Amplitudes larger 
than Px = 15 G cannot be used because the system col¬ 
lapses, leading to contacts between the beads. Close to 
collapse (regime IH), speed exhibits large fluctuations, 
as seen in Figure l^a). It has proven hard to control 
the speed and the pulsation mode in that regime. We 
also noticed that large offset values reduce the efficient 
locomotion (regime II) into a sharp range of /3x values. 

In order to quantify the triplet deformations, we com¬ 
puted the standard deviation 6 of the internal angle de¬ 
partures {ai — 60) from the regular triangle (in degrees), 
i.e. 
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The maximum value of this parameter over a period 
max {(5} is shown in Figure [^b) as a function of /3x. 
When Px increases from zero, max {(5} presents a jump 
similar to v above Bxo before saturating to around 2 de¬ 
grees (regime II). The speed v and deformation 5 are 
therefore correlated. It should be stressed that the defor¬ 
mation is relatively weak such that the locomotion should 
be linked to some additional motion of the beads such as 
the structure rotation. However, for higher amplitudes 
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FIG. 4. Full control of swimming trajectories - The 

remote control of a microswimmer, by adjusting the horizon¬ 
tal field orientation e, allows us to follow various paths such 
as U-turn, corner and loops. Laid next to one another, the 
three letters of “ULg” are then obtained. Bead centers are 
tracked for drawing the triangles at each period, emphasizing 
the orientation of the structure. The colors indicate the time 
evolution of each trajectory from dark to red. 


close to the collapse (regime III), the triangle is highly 
stretched providing much higher values of S where the 
speed fluctuates. 

The major ingredient triggering the locomotion is the 
field amplitude /3a,. Once the triangle pulsation is se¬ 
lected, the oscillation mode remains unchanged and the 
swimmer trajectory is a straight line. Although different 
from the x-axis, the swimming direction is in fact cor¬ 
related to the oscillating field direction. Indeed, Figure 

presents three independent trajectories of swimmers 
which are prepared by adjusting slowly the horizontal 
magnetic field direction in the x — y plane in order to 
induce U-turn, corner and loops, proving that the re¬ 
mote control can be fully exploited. The paths form the 
letters of “ULg”, being the logo of our University. The 
video of the U-turn is given in the supplementary ma¬ 
terials. Slowly adjusting the horizontal field orientation 
offers therefore a convenient way to control the path of a 
microswimmer. 

Discussion 

Two main mechanisms can provide the breaking of 
time reversibility necessary for low Reynolds locomotion. 
Symmetry breaking can either occur in the sequence of 
shapes adopted by the swimmer [1] or in the viscous 
drag felt by the swimmer m- Here, comparing typi¬ 
cal strength of magnetic force ^ and hydrody¬ 
namic interaction ~ ^ , where m is the magnetic 

moment of a bead, R its radius, U its speed and r is the 
center-to-center distance, yields a ratio ^ ~ 5 10^. The 
shape of the swimmer is thus driven by magnetic dipole- 
dipole interactions in our case. 

A magnetocapillary model, similar to m and based 


on the magnetic and capillary interactions only, allows 
us to capture the possible beads configurations when 
both vertical and horizontal magnetic fields are applied. 
Moreover, it provides clues for the origin of the non¬ 
reciprocal motion behind low Reynolds locomotion, as 
seen below. For convenience, distances rij between beads 
i and j are adimensionalized by the capillary length, 
being the c haracteristic length of capillary interactions 
^ = \/l/P9 ~ 2.7 mm [14]. Considering only quasi-static 
situations, i.e. excluding hydrodynamic interactions, the 
dimensionless bead-bead potential is given by 
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where Mcz and Mcx are magnetocapillary numbers [20| 
measuring the competition between magnetic and capil¬ 
lary interactions. The Bessel function Ko{r) is assumed 
to capture the attractive capillary interaction El- Each 
magnetocapillary number is proportional to dipole-dipole 
interactions. One can fix the number Mcz in order to ob¬ 
tain an equilibrium distance Vij = 2D/X = 0.3 when the 
horizontal field is zero (Mcx = 0, Mcz ^ 0.01). This 
distance corresponds to what we observe in our exper¬ 
iments. The parameter ^Mcx/Mcz becomes then the 
only relevant parameter in the problem, reducing to the 
ratio BxjBz. In order to study the equilibrium configu¬ 
rations, we start from a regular triangle with edge length 
r = 0.3, the horizontal field is slightly modified and the 
new equilibrium situation is found by simulated anneal¬ 
ing searching for the local minimum ofu = 1 ^ 12 +'^^ 23+'^315 
by producing tiny random moves for bead positions. The 
results of these simulations are discussed below. 

Figure [^a) presents the evolution of the internal an¬ 
gles ai of the triplet at equilibrium and for different val¬ 
ues of an increasing then decreasing horizontal field, i.e. 
when the ratio B^/B^ follows a quasi-static cycle. Ar¬ 
rows indicate how the structure evolves. The collapse 
of the structure takes place around B^jB^ ~ 1/2, as 
expected. In presence of a horizontal field, the trian¬ 
gle is roughly isosceles during the whole cycle since two 
angles are always close together. At low and increas¬ 
ing field strengths, one angle is larger than 60° and the 
other angles have similar values below 60°. For the lack 
of a proper term, we labelled this first type of isosceles 
“platy-isosceles”. It is characterized by internal angles 

} ( 3 ) 

and occurs when one edge is perpendicular to the field 
axis, as sketched in the inset of Figure E a). By sym- 
metry, two configurations can be obtained : left (<l) and 
right (>). Above some field strength, the system drasti¬ 
cally changes since two angles have similar large values 
while the later is much smaller than 60°. One has a sec¬ 
ond type, labelled “lepto-isosceles”, given by 

, 60 - 72(5 
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FIG. 5. Non-reciprocal motion - (a) Internal angles given 
by the model as a function of the increasing or decreasing 
horizontal magnetic interaction (see arrows). Angles ai are 
colored after being sorted from the smallest to the largest. 
Two triangles are illustrated using the same color code in or¬ 
der to evidence the orientation of the structure at two specific 
situations : < and V states, (b) Three periods T of the an¬ 
gles ai using the same color code in the real experiments. 
One observes that the blue curve oscillates between the red 
and green ones in each period. The system switches peri¬ 
odically between < and V isosceles, (c) The angle 0 of the 
whole structure during three periods. Oscillations are seen 
evidencing the successive switches, (d) This sketch shows the 
isosceles states as a function of 0 forming a loop at each cycle, 
i.e. a non-reciprocal motion. 


which is encountered when one side is parallel to the 
field. Two equivalent configurations are found : up (a) 
and down (v). In fact, this bifurcation between platy to 
lepto-isosceles, e.g. <1^ V, is accompanied by a rotation 
of 30° of the entire structure, as sketched in the Figure 
l^a). The system remains trapped in this last configura¬ 
tion till returning to a regular triangle at zero field. Such 


a switching behavior between two states is observed peri¬ 
odically in experiments. Indeed, Figure [^b) presents the 
evolution of the angles ai during three successive periods 
using the same color code : largest angle in red, lowest 
angle in green. During each period, the main qualita¬ 
tive characteristics predicted by the model are recovered 
: short successive periods during which both types of 
isosceles are seen. It should be noted that the model does 
not predict which configuration is chosen when the field 
goes back to zero. At zero field, small residual magneti¬ 
zation becomes relevant and might explain why the same 
configuration is always chosen, such that deformation is 
periodical. On should also remark that the experimental 
curves are much more complex than numerical ones due 
to the presence of hydrodynamic interactions and addi¬ 
tional effects, such as a slow rotation of the beads due to 
tiny residual magnetization. Because of the presence of 
these additional effects, the angles exhibit periodic jumps 
from one isosceles type to the other one. Figure [^c) 
presents the evolution of the orientation of the entire 
structure during the same three periods. The angle 0 
is seen to oscillate over one period, as expected. Figure 
id) illustrates the non-reciprocity of the shape succes¬ 
sion necessary to propulsion. It has a clear magnetocap¬ 
illary origin since each magnetic cycle induces a loop in 
the space of configurations. The successive switches be¬ 
tween isosceles states are indeed the driving mechanism 
of locomotion. 

Mesoscale self-assembly [9] is recognized as an elegant 
way to fabricate microsystems. Because no microfabri¬ 
cation techniques are involved, the experimental setup is 
quite straightforward. Furthermore, the use of magnetic 
fields to control the assembly and power the motion offers 
great flexibility. A magnetocapillary swimmer can indeed 
be assembled, adjusted, destroyed and reassembled as of¬ 
ten as necessary by only changing the magnetic fields. 
The spontaneous organization of a remote-controlled mi¬ 
croswimmer thus represents a valuable achievement. Al¬ 
though the motion of the self-assembled swimmer is con¬ 
strained along a liquid interface, many applications can 
be cited : for instance cargo transport, fluid mixing or mi¬ 
cromanipulator. Complex functionalities can be reached 
by considering specific beads and/or a larger number of 
beads. The fact that the setup does not require the use 
of chemically active substances is an advantage in this 
regard. Besides applications, this bead triplet has the 
great advantage to possess a low number of degrees of 
freedom such that fundamental aspects of low Reynolds 
locomotion ED [22] can now be experimentally explored 
in various conditions. 

Methods 

Ffereinafter, we present methods and experimental setup. 
A large Petri dish is filled with water. The liquid/air in¬ 
terface is placed at the center of a triaxial Helmholtz 
system that compensates the Earth magnetic field, and 
that is able to produce a uniform field in any direction. 
When a current i is injected in such coils, a uniform and 
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vertical magnetic field Bz is obtained in the Petri dish. 
Magnetic fields up to 30 G have been considered. Oscilla¬ 
tions of the horizontal field are provided by a function 
generator and an amplifier. Chrome steel particles (se¬ 
lected alloy AISI 52100, ps = 7830 kg/m^) are soft ferro¬ 
magnetic beads, and they do not exhibit any hysteretic 
behavior in the range of field values used herein. As 
a result, particles retain a negligible residual magnetic 
moment once the field is removed. We have estimated 
this residual magnetization to be roughly 50 times lower 
than the induced magnetization due to the applied field. 
Prior to experiments, spheres are washed with isopropyl 
alcohol and thereafter dried in an oven. Different bead 
diameters have been studied but the results shown herein 
correspond to D = 500 pvn. At this scale, partial wetting 
ensures the floatation of the spheres. A CCD camera 
records images from above. Image analysis provides the 
position of each bead as a function of time. 
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